OPEN 3 ACCESS Freely available online 



•0-PLOS I ONE 



Distribution and Frequency of kdr Mutations within rm\ 
Anopheles gambiae s.l. Populations and First Report of 
the Ace. 7G119S Mutation in Anopheles arabiensis from 
Burkina Faso (West Africa) 

Roch K. Dabire 1 *, Moussa Namountougou 1 , Abdoulaye Diabate 1 , Dieudonne D. Soma 1 , Joseph Bado 1 , 
Hyacinthe K. Toe 1,2 , Chris Bass 3 , Patrice Combary 4 

1 IRSS (Institut de Recherche en Sciences de la Sante), Centre Muraz, Bobo-Dioulasso, Burkina Faso, 2 Department of Vector Biology, Liverpool School of Tropical Medicine, 
Liverpool, United Kingdom, 3 Biological Chemistry and Crop Protection Rothamsted Research, Harpenden, United Kingdom, 4 National Malaria Control Programme, 
Ministry of Health, Ouagadougou, Burkina Faso 



Abstract 

An entomological survey was carried out at 15 sites dispersed throughout the three eco-climatic regions of Burkina Faso 
(West Africa) in order to assess the current distribution and frequency of mutations that confer resistance to insecticides in 
An. gambiae s.l. populations in the country. Both knockdown (kdr) resistance mutation variants (L1014F and L1014S), that 
confer resistance to pyrethroid insecticides, were identified concomitant with the ace-1 G119S mutation confirming the 
presence of multiple resistance mechanisms in the An. gambiae complex in Burkina Faso. Compared to the last survey, the 
frequency of the L1014F kdr mutation appears to have remained largely stable and relatively high in all species. In contrast, 
the distribution and frequency of the L1014S mutation has increased significantly in An. gambiae s.l. across much of the 
country. Furthermore we report, for the first time, the identification of the ace.l G116S mutation in An. arabiensis 
populations collected at 8 sites. This mutation, which confers resistance to organophosphate and carbamate insecticides, 
has been reported previously only in the An. gambiae S and M molecular forms. This finding is significant as 
organophosphates and carbamates are used in indoor residual sprays (IRS) to control malaria vectors as complementary 
strategies to the use of pyrethroid impregnated bednets. The occurrence of the three target-site resistance mutations in 
both An. gambiae molecular forms and now An. arabiensis has significant implications for the control of malaria vector 
populations in Burkina Faso and for resistance management strategies based on the rotation of insecticides with different 
modes of action. 
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Introduction 

The pyrethroid class of insecticides have become a mainstay for 
vector control since the ban of DDT due to off-target toxicity and 
the development of resistance. They have been most widely used 
to treat bed nets (ITNs) dedicated to personal and community 
protection [1,2,3]. Unfortunately, knock down resistance (kdr) to 
pyrethroids, which also confers cross-resistance to DDT, was first 
reported in Anopheles gambiae populations from Cote d'lvoire [4]. 
Resistance likely resulted from the earlier intensive use of DDT 
and selection from pyrethroid use in crop protection particularly in 
cotton areas [5,6]. kdr was initially shown to result from a point 
mutation (L1014F) in the pyrethroid target protein the voltage- 
gated sodium channel [7]. Based on a simple PCR diagnostic 
developed in the first report of the kdr mutation [7] several studies 
have been carried out on the distribution and the frequency of this 



mechanism throughout Africa. Initial studies showed that L1014F 
kdr was most widely distributed in West African An. gambiae s.l. 
populations [6,8,9]. This mutation was observed initially in the S 
molecular form of An. gambiae s.s. reaching high frequency but 
was not found either in sympatric mosquitoes of the M molecular 
form or An. arabiensis populations [5]. This provided further 
evidence of reproductive barrier between the M and S molecular 
forms [10,11] and the two molecular forms of An. gambiae s.s. 
were recendy confirmed as two distinct species termed Anopheles 
coluzzii for the M form and Anopheles gambiae for the S form [12]. 
However, a few years after the initial finding of the kdr mutation in 
the S molecular form, this mutation was also reported in the M 
form from the littoral of Benin and Cote d'lvoire [13]. In-depth 
investigations carried out later in these geographic regions 
confirmed that this phenomenon was frequently observed in 
littoral but was rare inland [11]- DNA sequencing of these 
mosquitoes suggested that the mutation emerged in the M form by 
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genetic introgression from the S form [14,15]. In contrast, the 
emergence of the Leu-Phe kdr mutation within Anopheles 
arabiensis resulted from a de novo mutation event [15]. An 
extensive monitoring program in Burkina Faso has revealed that 
the L1014F kdr mutation initially detected in low frequency in the 
An. gambiae M molecular form and An. arabiensis [11,15] has 
spread throughout the country and is observed in mosquito 
populations at relatively high frequency [16,17], Recently the 
L1014S kdr, which initially predominated in East Africa [18,19], 
was reported in West Africa, first in Benin and then Burkina Faso 
within An. arabiensis populations [20,21]. More recendy this 
mutation was reported in a small number of individuals of the M 
and S forms otAn. gambiae in Burkina Faso [22]. Taken together 
these results provide fundamental insight into the evolutionary 
processes underlying resistance in Anopheles gambiae s.l. Further- 
more from an applied perspective, the emergence of resistance has 
significant implications for vector control programmes, especially 
those focused on the use of ITNs/Long-Lasting Insecticidal Nets 
(LLINs) or indoor residual sprayings (IRS). Although LLINs had 
shown good control of certain pyrethroid resistant populations 
[23] reduced efficacy of treated nets against An. gambiae 
populations with kdr resistance has since been reported [24]. 

Other insecticides belonging to the organophosphate (OP) and 
carbamate (CM) classes have been investigated to be used in 
mosaic, or in combination, with pyrethroids for bednet impreg- 
nation [25]. In addition to the use of LLINs, bendiocarb was 
recently used in IRS applications in West Africa through the 
President's Malaria Initiative (PMI) roadmap [26]. Initially 
described in Culex populations from Cote-d'Ivoire [27] reduced 
susceptibility to OPs and CMs was observed in An. gambiae 
populations in the North of Cote d'lvoire and related to the 
domestic use of insecticide [28]. An. gambiae populations from 
Benin with resistance to the CM bendiocarb were reported after 
just three year of IRS use [29]. A common mechanism of 
resistance to OP and CM insecticides results from a single point 
mutation (termed ace-l R )in the target protein the acetylcholines- 
terase enzyme [30]. This mutation results in a glycine to serine 
replacement at amino acid position 1 19 and can be detected by a 
simple PCR-Restriction Fragment Length Polymorphism (RFLP) 
diagnostic [31]. This approach has been used to examine the 
frequency and distribution of this mutation in Burkina Faso where 
it was found predominately in the An. gambiae S form and in low 
frequency in the M form [9,16,32]. A recent study suggested that 
the mutation had introgressed from one form to the other but the 
precise origin of the introgression could not be determined due to 
the small sample size [33]. Since then, extensive country- wide 
surveys were performed in Burkina Faso from 2008 to 2010 and 
no case of An. arabiensis carrying this mutation was reported, 
although sample sizes for this species were sometimes small 
[16,17]. 

However insecticide resistance may also occur by other 
physiological mechanisms such as metabolic detoxification 
through increased enzyme activities (monooxygenases, esterases, 
or glutathione S- transferases) [34,35]. 

Burkina Faso is composed of three agro-climatic areas which 
exhibit different patterns of insecticide use especially in relation to 
crop protection. The present study provides an update on the 
distribution and the prevalence of the kdr L1014 and L1014S 
andace-l R mutations in An. gambiae s.l. populations throughout 
the 1 3 health regions dispersed across these different agro-climatic 
areas. We report here, for the first time, the occurrence of the ace- 
1 R mutation at remarkably high frequencies in An. arabiensis. 



Materials and Methods 

Study sites 

Burkina Faso covers three ecological zones, the Sudan 
savannah zone in the south and west where rainfall is relatively 
heaviest (5-6 months), the arid savannah zone (Sudan-sahelian) 
which extends throughout much of the central part of the 
country and the aridland (Sahel) in the north. The northern part 
of the country has a dry season of 6—8 months. The varied 
ecological conditions are reflected in the different agricultural 
systems practiced throughout the country, from arable to pastoral 
lands. The western region constitutes the main cotton belt 
extending to the south where some new cotton areas have been 
cultivated since 1996. All ecological zones support the existence 
of Anopheles species that vector malaria and the disease is 
widespread throughout the country. Larvae were sampled from 
15 sites dispersed throughout the three ecological zones 
(Table 1). The GPS coordinates were incorporated in Table 1. 

Mosquito sampling 

Larvae of An. gambiae s.l. were collected from at least 10 
breeding sites dispersed throughout each sampling site mainly 
comprising pools of standing water and other small water 
collections. Larvae were pooled to constitute a colony, which was 
reared in the insectary to adulthood. A sample of 100 adult 
females were randomly sorted, killed and kept on silica gel in 1.5- 
ml tubes and stored at — 20°C prior to PCR analysis. Anopheline 
species were identified morphologically using the standard 
identification keys of Gillies and Cootzee [36] . 

PCR analyses 

An average of 30 mosquitoes was sampled per site by PCR 
analysis. Genomic DNA was extracted from single specimens and 
used as template for PCR to determine the species within the An. 
gambiae complex using the protocol SINE 200 of Santalomazza et 
al. [37] that allows the concomitant identification of An. gambiae 
M and S (respectively known as Anopheles coluzzii and Anopheles 
gambiae) and An. arabiensis. The same individuals were then 
tested for both the L1014F and L1014S kdr mutations using the 
protocols of Martinez-Torres et al. [7] (using specific primers Agdl, 
Agd2, Agd3 and Agd4) and Ranson et al. [18] (using Agdl, Agd2, 
Agd4 and Agd5) respectively: 

- Agdl: 5'-ATAGATTCCCCGACCATG-3'; 

- Agd2: 5 ' -AGAC AAGGATGATGAAC C-3 ' ; 

- Agd3: 5'-AATTTGCATTACTTACGACA-3'; 

- Agd4: 5 ' -CTGTAGTGATAGGAAATTTA-5 ' ; 

- Agd5: 5 '-TTTGC ATTACTTACGACTG-3 ' . 

The ace-l R mutation was detected from the same samples by 
PCR according to the protocol of Weill et al. [31] using specific 
primers Ex3AGdir (GATC GTGG AC AC C GTGTTC G) and Ex3- 
AGrev (AGGATGGCCCGCTGGAACAG). Then the PCR 
products were digested using Alu 1 enzyme at 37°C for 3 hours. 

Statistical analysis 

Data were compared between ecological zones and pooled for 
each species to compare the genotypes frequency between An. 
gambiae species by Chi 2 tests. The genotypic frequencies of 
L1014F and L1014S and ace-l R in mosquito populations were 
compared to Hardy-Weinberg expectations using the exact test 
procedures implemented in GenePOP (ver.3.4) software [38]. 
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An. An. k Sudan Sudan- Sahel 



gamhiae ro/urrn arabiensis sahcl 

Figure 1. Comparison of allele frequencies of 1014F, 1014S and ace-f mutations within Anopheles gambiae, An. coluzzii and /4n. 
arabiensis populations from 1 5 sites dispersed across the 3 agro-ecological regions of Burkina Faso. 

doi:10.1371/journal.pone.0101484.g001 



Ethical issues 

Ethical approval was not required in this study. 

This study was not carried out on private land. For each, no 
permission was required our study does not degrade the 
environment. No permission was required for these locations/ 
activities as the field activities did not involve damaged of 
protected species. We did not use any vertebrate during this study. 

Results 

Out of 5 1 6 mosquitoes analysed in PGR, 513 successfully scored 
(less than 5% failure rate). Overall species composition of the 
collected mosquitoes comprised a higher proportion of An. 
gambiae (51.7%) than An. coluzzii (21.6%) and An. arabiensis 
(26.7%) (Table 1). The species repartition across the three 
ecological regions revealed that An. gambiae was the predominant 
species in all regions including, in the Sahel where it comprised 
more than 49% of the An. gambiae s.l. population. Anopheles 



arabiensis was the second most predominant vector found in 
samples collected from the three regions. Somewhat An. coluzzii 
was found at a relatively low proportion of less than 15%. The 
central areas were characterised by an overlapped repartition of 
the three species 38.4%, 27.81% and 33.75% for An. gambiae, 
An. coluzzii and An. arabiensis respectively and proportions did 
not differ significantly [% 2 = 1.95, df= 1, P>0.05). In the Sahel 
region, An. gambiae also predominated (49.75%) and the 
proportions of the two other species did not differ significandy at 
21.01% and 29.74% for An. coluzzii and An. arabiensis 
respectively (% 2 = 4.88, df= 1, P>0.05). 

The overall frequency of the L1014F mutation averaged 50% 
and did not significantly differ between species (Figure 1A) 
whatever the ecological zone (Figure IB) (x 2 = 0.14, df= 1, P> 
0.05) even though the highest values were observed in the sudan 
zone (Figure 2). However some deviation from Hardy- Weinberg 
expectations was observed within the An. arabiensis populations in 
Dedougou and Dori and within An. coluzzii populations in Fada, 
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Figure 2. Distribution the 1014F kdr allele frequency from 15 sites dispersed across Burkina Faso. 

doi:1 0.1 371 /journal.pone.01 01 484.g002 



Kaya, Ouahigouya and Dori with an excess of resistant 
homozygous alleles (Table 2). The same patterns were found in 
seven sites for An. gambiae (Gaoua, Banfora, Sindou in the West, 
Dedougou, Koudougou and Koupela in the central region and 
Ouahigouya in the Sahel) (P<0.05). 

The overall allele frequency of the L1014S kdr mutation 
(Figure 3) was relatively higher in An. gambiae (48%) followed by 
An. coluzzii (38%) and An. arabiensis populations (37%) with no 
significant difference between the last two (% 2 = 3.24, df= 1, P> 
0.05) (Figure 1C). Comparing between ecological regions, L1014S 
kdr frequency did not differ significantly between species, except in 
the Sahel where it was significantly higher in An. coluzzii than An. 
arabiensis (x 2 = 10.21, df= 1, P<0.001) and An. gambiae (P< 
0.04) (Figure ID). The observed genotypic frequencies were not 
significantly different from Hardy-Weinberg expectations at the 
95% confidence level (Table 2) in populations from any site except 
in the An. gambiae populations from Orodara, Soumousso, 
Koupela, Fada, and in the An. arabiensis populations from 
Dioulassoba and Kaya where a heterozygous deficit was observed 
(P = 0.005) and An. gambiae populations in two sites (Dedougou 
and Kaya) where an excess of heterozygotes was observed (P< 
0.05). 

The ace-l R mutation (Figure 4) was recorded in all the 15 sites 
under study with a wider distribution within the An. gambiae 
populations (Table 3). The overall allele frequency of ace-l R was 
significantly higher in An. arabiensis (0.26) than in An. gambiae 
(0.11) (x 2 =14.4; df=l, P = 0.001) and An. coluzzii (0.09) 



(X 2 = 11.77, df=l, P = 0.006) (Figure IE) with no significant 
difference between the last two (/ 2 = 0.37, df= 1, P = 0.54). 
Compared between zones, the ace-l R allele frequency in An. 
arabiensis was higher than that of An. coluzzii (x 2 = 8.15, df= 1, 
P = 0.004) and An. gambiae (x 2 = 9.79, df= 1, P<0.001) in the 
Sudan and Sudan-sahelian savannah (with respectively % 2 = 6.89, 
df= 1, P<0.008 and % 2 = 17.34, df= 1, P<0.0003) (Fig. IF). In 
the Sahel no significant difference was observed between the three 
species (x 2 = 0.89-0.021, df= 1, P>0.05). The observed genotypic 
frequencies were significantly different from Hardy-Weinberg 
expectations at the 95% confidence level (Table 3) in An. gambiae 
population from Orodara, Soumousso, Koudougou, Fada, 
Ouahigouya, Dori and Dioulassoba, Koudougou and Kaya for 
An. arabiensis where a heterozygote deficit was observed 
(P = 0.005). Furthermore, the percentage of homozygous resistant 
individuals was significantly higher in An. arabiensis (25%) than in 
An. gambiae (6.25%). No homozygous resistant individual was 
recorded in An. coluzzii from any site. 

Discussion 

This study provides current information on the distribution of 
three members of the Anopheles gambiae complex across Benin 
and the frequency and distribution of three important target-site 
resistance mechanisms in these populations. In regards to the 
distribution of An. gambiae species throughout the country, the 
most significant finding is that An. arabiensis appears to be 
spreading in the Sudan whereas in the past it comprised only 
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Figure 3. Distribution the 1014S kdr allele frequency from 15 sites dispersed across Burkina Faso. 

doi:1 0.1 371 /journal.pone.01 01 484.g003 



around 5% of the An. gambiae complex species [6]. Furthermore, 
this species is now present in Sindou at 14.29% (nearest the 
frontier of Cote-d'Ivoire) where it was absent a decade ago [9] . 
The reason for this is not clear but could be related to climatic 
changes, such as irregularities in rainfall observed in the 
boundaries of the Sudan region that may make the landscape 
more favourable to the establishment of this species. 

Across sampling covering 1 5 sites we identified the L 1 0 1 4F and 
L1014S Mr mutations concomitant with the ace- 1 G119S 
mutation confirming the presence of multiple resistance mecha- 
nisms in the An. gambiae complex in Burkina Faso [16,17]. The 
distribution and the prevalence of the L1014F kdr mutation in An. 
gambiae species including An. gambiae, An. coluzzii and An. 
arabiensis, has been well documented in Burkina Faso for over a 
decade [9,16]. Many studies reported this mutation at high 
frequency within An. gambiae and An. coluzzii populations 
especially in An. gambiae populations from the Sudan area where 
mutation frequency was approaching fixation [9,15,16]. Over 
recent years the frequency of this mutation has increased within 
both An. coluzzii and An. arabiensis. In this study although the 
L1014F mutation remains widespread in all three ecological 
regions and is present at relatively high frequency within the three 
species (averaging 50%), the frequencies reported in this current 
study were lower in the Sudan ecological regions (West and South 
West covering the old cotton belt) than those from previous studies 
[9,16,22]. For the other climatic zones i.e. central and northern 
regions the allele frequencies of L 1 0 1 4F varied within the three 



species with particularly high frequencies in An. arabiensis. The 
reason(s) for the reduction of L1014F frequency in An. gambiae 
populations in the Sudan area is not known, however, a similar 
trend was recently observed in the Western region of Burkina Faso 
where transgenic and biological control practices have been 
implemented for crop protection of cotton over the last four years 
(a long side conventional crop protection approaches) (Namoun- 
tougou, unpublished). These alternative cotton-growing practices 
would be expected to reduce the quantity and frequency of 
insecticide use in agriculture and this may in turn reduce the 
selection pressure experienced by local mosquito populations. The 
analysis of observed genotypic frequencies revealed a heterozygote 
deficit for the L1014F mutation in the three species of^4n. gambiae 
s.l. from many sites especially in the Sahel for An. coluzzii and An. 
arabiensis and in the Sudan and Sudan-Sahel for An. gambiae 
which deviated significantly from Hardy-Weinberg expectations. 
This finding is not surprising as the same patterns were observed in 
the West (Orodara and Soumousso) four years ago [9] in 
combination with a novel mutation, N1575Y, in the voltage-gated 
sodium channel, recently reported in An. gambiae s.l. populations 
in Soumousso [39]. 

The L1014S kdr mutation was recently recorded at highest 
frequency in An. arabiensis populations in the centre on the 
country [21] and in Bobo-Dioulasso at frequencies averaging 38% 
[40]. Previous studies have recorded only a few individuals of An. 
gambiae and An. coluzzii from the Centre-East part of the country 
[17] carrying this mutation in the heterozygous form. The present 
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Figure 4. Distribution the ace- 1 R allele frequency from 15 sites dispersed across Burkina Faso. 

doi:10.1371/journal.pone.0101484.g004 



study reveals that this mutation has since spread across the whole 
country and is now observed at relatively high and similar 
frequencies (40%) between the three species. The comparison of 
the observed genotypic frequencies of this mutation with that 
expected for Hardy-Weinberg equilibrium indicated, depending 
on the site, a deficit or excess of heterozygotes, mainly for An. 
gambiae populations. The occurrence of the L1014F kdr mutation 
in An. coluzzii had been suggested to have occurred by 
introgression from An. gambiae and via a de novo mutation event 
in An. arabiensis [15], however, the origin of the L1014S mutation 
in An. gambiae, An. coluzzii and An. arabiensis species in West 
Africa is not so clearly understood. The proximity of Burkina Faso 
from the Benin frontier where the L1014S mutation was first 
reported in An. arabiensis populations [20] suggests that it arrived 
in Burkina Faso via migration of An. arabiensis carrying the 
mutation from Benin, however, the origin of this mutation in An. 
gambiae and An. coluzzii populations in Burkina Faso remains to 
be elucidated. 

In this study we report, for the first time, the presence of the 
ace.l G119S mutation in An. arabiensis populations from eight 
sites: Dioulassoba, Soumousso in the West, Boromo, Dedougou, 
Koudougou, Nanoro and Fada in the Centre-North and East and 
Kaya in the North. In these sites An. arabiensis was observed as 
the second major vector after An. gambiae except at Fada and 
Nanoro where the proportion of An. arabiensis was lower than 
that of An. coluzzii. To confirm this finding, we repeated the PGR 
amplification of ace. / R for our An. arabiensis specimens and used, 



as a control, 30 specimens of An. Arabiensis which we had 
confirmed in a previous study do not have this mutation. No false 
positives were observed in these samples suggesting our data is 
robust. The ace. 1 allele was observed in this study in An. 
arabiensis at varying frequency reaching a maximum value of 78% 
in populations from Dioulassoba and the lowest value in Kaya at 
8%. Except for samples from Soumousso and Nanoro where the 
sample size was not sufficient (n<10) to compare genotype 
frequencies, deviations from Hardy-Weinberg equilibrium were 
observed at three sites (Dioulassoba, Koudougou and Kaya) as a 
result of a high heterozygote deficit. The same pattern was 
observed in An. gambiae from Orodara, Soumousso, Koudougou, 
Fada, Ouahigouya and Dori. The deficit of heterozygous 
genotypes observed in Orodara and Soumousso is not new as 
Dabire et al. [41] reported similar results from the these areas from 
which the duplicated allele (ace.l D ) was reported by Djogbenou et 
al. [33]. It is possible that this duplicated allele ace.l D is also 
present within An. arabiensis especially in Dioulassoba where the 
proportion of homozygous mutants was atypically high (60%). The 
high frequency of this mutation in Dioulassoba populations is 
intriguing as recent studies failed to find any L 1 0 1 4F kdr or ace-l R 
in An. arabiensis population from this site [40,42]. As for the 
L1014S mutation, additional sequence analysis of the region 
flanking the ace.l locus are necessary to confirm whether the ace.l 
mutation in An. arabiensis has evolved along the same pathway as 
kdr e.g. as a de novo mutation or introgression from An. gambiae 
or An. coluzzii. Unfortunately our PCR data is not backed up by 
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insecticide susceptibility bioassays and so we cannot assess the 
correlations between kdr and ace- 1 mutations and the phenotypic 
expression of resistance. 

The emergence of the ace-l R mutation in An. gambiae s.l. 
population from the cotton-growing areas may be linked to the 
agricultural use of OP and CM insecticides used for crop 
protection. Other sources of selection pressure outside the cotton 
belt include insecticide use for vegetable growing and domestic use 
of insecticide in public health. Bioassays performed in 2012 on An. 
gambiae populations from sites located in the cotton belt of the 
West of Burkina Faso revealed the development of resistance to 
CMs and OPs especially to benidocarb (Dabire, unpublished) 
correlating with the prevalence and frequency of genetic resistance 
revealed in the present study. However, further bioassays on a 
wider scale are now required in order to understand the 
implications of the current status of the ace-l R mutation for the 
efficacy of OP and CM insecticides in vector control in Burkina 
Faso. The information provided by such studies combined with 
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